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Abstract 

We study the relation between the thermodynamics and held equations of generalized gravity 
theories on the dynamical trapping horizon with sphere symmetry. We assume the entropy of 
dynamical horizon as the Noether charge associated with the Kodama vector and point out that it 
satisfies the second law when a Gibbs equation holds. We generalize two kinds of Gibbs equations 
to Gauss-Bonnet gravity on any trapping horizon. Based on the quasi-local gravitational energy 
found recently for f(R) gravity and scalar-tensor gravity in some special cases, we also build up the 
Gibbs equations, where the nonequilibrium entropy production, which is usually invoked to balance 
the energy conservation, is just absorbed into the modified Wald entropy in the FRW spacetime 
with slowly varying horizon. Moreover, the equilibrium thermodynamic identity remains valid for 
f(R) gravity in a static spacetime. Our work provides an alternative treatment to reinterpret the 
nonequilibrium correction and supports the idea that the horizon thermodynamics is universal for 
generalized gravity theories. 
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I. INTRODUCTION 



With the discovery of quantum Hawking radiation, it became clear that a black hole 
behaves as an ordinary thermodynamic system with temperature proportional to surface 
gravity and entropy measured by its horizon area. The laws of black hole mechanics which 
are implied in Einstein field equations can be treated as the laws of thermodynamics. To 
comprehend why gravity knows thermodynamics, Jacobson turned the logic around and 
disclosed that Einstein's equation can emerge as an equation of state from the basic ther- 
modynamic relation in local Rindler spacetime This puzzling thermodynamic feature 
of gravity and/or spacetime acted as an important motivation of the proposal that grav- 
ity might not be a fundamental interaction but rather an emergent large scale/numbers 
phenomenon 2|, y] . If this were true, not only general relativity but also more generalized 
theories of gravity, such as the ones with higher order curvature corrections, should be seen 
under this same light, presuming higher order curvature terms originated from the quantum 
corrections of underlying microscopic theory. 

Actually, it was disclosed by Padmanabhan et al. that beside Einstein's equation 4, 
[5(1, the field equation of Gauss-Bonnet gravity (as well as for a wider class of Lanczos- 
Lovelock gravity) can be written as an equilibrium thermodynamic identity [6], which is a 
Gibbs (-like) equation (or the so called first law of thermodynamics), near the horizon of 
a static spherically symmetric spacetime. This result has been extended to general static 
spacetimes recently 2fl. But for the dynamical spacetime, Padmanabhan's identity has 



not been genera 
unified first law 
FRW spacetime 



ized. There exists another similar Gibbs equation based on Hayward's 
8|, which has been constructed_on two special dynamical spacetimes, the 

Especially in the FRW spacetime, more 



10( and Vaidya spacetime [11 



12 



14] and loop quantum gravity [15j have 



modified gravity theories like braneworld gravity 
been similarly described. 

However, the equilibrium thermodynamics, which is suitable for the Gauss-Bonnet gravity 
and the mentioned other gravity theories, has not been extended trivially to f(R) gravity 
and scalar-tensor gravity. Instead, it was shown that a nonequilibrium entropy production 
term needs to be invoked to balance the energy conservation 16|, [l7|. In the local Rindler 



spacetime, the entropy production term is presumed originated from the bulk viscosity of 
internal spacetime. In this approach, it has been shown that even for Einstein gravity 
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an entropy production term may still appear, which is interpreted as the space-time shear 
viscosity with the ratio to entropy density as 1/ (47r), consistent with the result of AdS/CFT 
duality. In a dynamical FRW spacetime, even for the case with slowly varying horizon, there 
are also some extra terms when one tries to reconstruct the field equation as an equilibrium 
identity of thermodynamics {9], [l8^ 20 1 . These terms have been assumed as the entropy 
production as well but it is not clear whether it can be attributed to the spacetime viscosity 
too. 

One alternative treatment to reinterpret the nonequilibrium correction was studied in the 
FRW spacetime with slowly varying horizon 21], where a mass-like function is introduced 
to absorb the entropy production terms. This mass-like function in Einstein gravity has a 
close relation to Misner-Sharp energy (mass) 22], which is a widely accepted quasi-local 
gravitational energy in the spherically symmetric spacetime and takes role as the gravi- 
tational energy in the mentioned two Gibbs equations. We have showed [19] the similar 
relation existed between the mass-like function and the generalized Misner-Sharp energy in 
Gauss-Bonnet gravity 23], but it is not known whether it actually has the significance of 
gravitational energy, especially in f(R) gravity. 

Another alternative method was proposed by Elizalde and Silv a |24lj . they noticed that 
the entropy of stationary horizon is well defined by Wald entropy 25], which is a Noether 
charge associated with the horizon Killing vector, but it is less understood for the horizon 
entropy in a dynamical spacetime, where the Killing vector can not be found in general. 
Iyer and Wald proposed that one can approximate the metric by its boost-invariant part to 
"create a new spacetime" where there is a Killing vector. However, the obtained dynamical 
entropy is not invariant under field redefinition in general 26] . Elizalde and Silva showed that 
the equilibrium thermodynamics can derive field equations of f(R) gravity, provided that 
the dynamical horizon entropy is still of same form as the stationary case but the entropy 
variation is evaluated by its boost-invariant part at leading order. We have proved that 
the same method can be suitable for the scalar-tensor gravity j^. Although this method 
has not invoked the dynamical entropy with complete metric, it suggests that the correct 
dynamical entropy may be essential to understand the nonequilibrium entropy production. 

Recently, one important step to realize the equilibrium or nonequilibrium spacetime has 
been taken, which shows that the field equations for arbitrary diffeomorphism-invariant 
gravity theories can be obtained as an equilibrium state equation of Rindler horizon ther- 
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modynamics 2814301] . Moreover, it proves that the Wald entropy satisfies the second law 



of thermodynamics when the null energy condition is met. If this approach is perfect, it 
completed the implementation of Jacobson's proposal to express Einstein's equation as a 
thermodynamic equation of state, and the nonequilibrium entropy production is not nec- 
essary. However, since the approximate Killing field is invoked and the horizon entropy is 
assumed as Wald entropy, this result may be not applicable to dynamical spacetimes. 

Thus it is urgent to study the dynamical horizon entropy and the corresponding horizon 
thermodynamics. Besides constructing the entropy by a created Killing vector with the 
boost-invariant metric, it is natural to consider a Noether charge construction associated 
with certain special dynamical vector field. Unfortunately, if the vector is not a Killing 
vector, the construction is not unique [26], because an arbitrary exact form can be added 
in Lagrangian which affects the Noether current but not the field equation. Moreover, the 
Noether current and the Noether potential have the freedom to the addition of a closed 
form. In the absence of a deeper understanding of black hole entropy, it seems impossible 
to impose a fundamental criterion in defining the entropy of a nonstationary black hole. 



However, there is an important assistant criterion given in [31] that if an entropy expression 
that satisfies the second law can be defined, that would be a preferred definition. Of course, 
the dynamical entropy should recover the Wald entropy in stationary cases, and can be 
expected with a covariant form. 

This kind of entropy expression has been presented for Einstein gravity in spherical sym- 
metry spacetimes. It is well known that in these spacetimes, there is a preferred time 
direction given by the Kodama vector 32|, which is a natural dynamical analogue of a 
stationary Killing vector. Hayward hence proposed that the Wald entropy can be alterna- 



tively associated with Kodama vector 



33 



34j. For Einstein gravity, the dynamical horizon 



entropy, which has been called as Wald-Kodama entropy, has the same simple form Aj '4 
(A is the area of horizon surface) as for stationary black holes. This was also justified by 
evaluating the surface terms in a dual-null form of the reduced action in two dimensions 
351 ] . Importantly, it was proved that the surface entropy satisfies the second law in the 



dynamical spacetime 



33|. 



For modified gravity theory, however, the Wald-Kodama entropy has not been obtained. 
In this paper, we will give the general expression of Wald-Kodama entropy and evaluate it for 
Gauss-Bonnet gravity, f(R) gravity and scalar-tensor gravity. One can find that the former 
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has the same form as the stationary case, while for later two cases the forms are different. We 
will show that if Hayward's equilibrium thermodynamic identity can be built up using these 
entropy expressions, they can satisfy the second law under the null energy condition. We will 
try to construct the equilibrium thermodynamics both in Padmanabhan's and Hayward's 
approaches. For Gauss-Bonnet gravity, our construction will invoke the generalized Misner- 
Sharp energy which has been obtained in any spherical symmetry spacetime 23(. For f(R) 



gravity and scalar-tensor gravity, a very recent work [36| has disclosed that the desired 
generalized Misner-Sharp energy can not be always found and written in an explicit quasi- 
local form. Fortunately, such a form exists for f(R) gravity and scalar-tensor gravity in 
an FRW universe and exists for f(R) gravity in the static spherically symmetric solutions 
with constant scalar curvature. We will show that for Gauss-Bonnet gravity, the equilibrium 
thermodynamics can be derived by Padmanabhan's and Hayward's methods, respectively. 
For f(R) gravity and scalar-tensor gravity in an FRW universe, it is interesting to see that 
the equilibrium thermodynamics also holds on slowly varying horizon. For the static case, 
Hayward's unified first law is trivial. But based on the generalized Misner-Sharp energy, we 
can set up Padmanabhan's identity for f(R) gravity too. 



II. ENERGY OF TRAPPING HORIZONS 



Suppose the n-dimensional spacetime (M n , g^) to be a warped product of an (n — 2)- 
dimensional spherical symmetry space (i^ n _ 2 , jij) and a two-dimensional orbit spacetime 
(M 2 , h ab ) under the isometries of (K n _ 2 , 7^). Namely, the line element is given by 

ds 2 = h a bdx a dx b + r 2 {x)^ijdy % dy\ 

where r is the areal radius for an in — 2)-sphere i^ n -2- It is useful to locally rewrite the line 
element in the double-null coordinates as 

ds 2 = -2e~^ v) dudv + r 2 {u, v)d 2 Q n _ 2 (1) 

where d 2 VL n _ 2 denotes the line element of the (n — 2)-sphere. The key geometrical objects 
are dynamical trapping horizons [3^], which are hypersurfaces H (in space-time) foliated 
by marginal surfaces. A marginal surface is a spatial surface on which one null expansion 
vanishes. The null expansions of two independent future-directed radial null geodesies are 
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expressed as 6 + = (n — 2)r~ 1 r,„ and 6L = (n — 2)r~ 1 r,t,. A marginal sphere with 9 + = is 
called future if 6L < 0, past if 6L > 0, bifurcating if 6L = 0, outer if d v 9 + < 0, inner if d v 6 + 
> and degenerate if d v 9 + = 0. We will follow Hayward's local definition of black holes (or 
white holes) as an outer trapping horizon and the black hole's possible inner boundary or 
cosmological horizon is taken as an inner trapping horizon. Recalling that in the spherically 
dynamical case, it is possible to introduce the conserved Kodama vector, given by 

K * = - e ^V,r = (e*d v r, -e+d u r, 0, 0, • • • ), (2) 

where e^ u = e a b {dx a ) (rfx 6 )^, e & is a volume element of (M 2 , h a b), and the minus in the 
definition is added to reduce the Kodama vector to the Killing vector = (1, 0, • ■ ■ ) for 
a static Schwarzschild spacetime. The dynamical surface gravity |8j associated with the 
trapping horizon can be defined directly from the Kodama vector j32| 

k = ~e ab \7 a K b = -e*d u d v r. (3) 

Similar to the Killing vector in the stationary case, the Kodama vector can be also related to 
certain gravitational energy and the entropy of trapping horizon. Let's consider the energy 
first in this section. One can define a current by the energy-momentum tensor of matter 
= —T^ U K V . If there is a conservation law V M = 0, an associated charge can be obtained 

as 

E = j J&IEp (4) 

by integrating the locally conserved currents over some spatial volume S with boundary. 

For Einstein gravity, the field equation is G^ v = 8ttT^ u , where we have set the gravita- 
tional constant G = 1. It is easy to prove that the current is conserved and the conserved 
charge is just the Misner-Sharp energy. 

For Gauss-Bonnet gravity with Lagrangian 

L = R + aL GB) (5) 

where a is the coupling constant and 

Lgb = R 2 — AR iW B! w + R flu \ p R' J ' uXp , 



the field equation is 



Gp, + aH^ = (6) 



where 



Hfj, u — 2(RR fJiU — 2R^\R^ — 2R Xp R fl \ up + R* pa R v \ pcT ) — -g pu L>GB- 



It has been proved that the current is still conserved, and the generalized Misner-Sharp 
energy has been obtained as [23] 

(n - 2)K- 2 r n - 3 



E 



16tt 



[(1 + 2e*r, u r, v ) + «r" 2 (l + 2e*r, u r, v ) 2 ] , 



(7) 



where V n -i denotes the area of (n — 2)-sphere K n _ 2 and a = (n — 3)(n — 4)<x 

For the nonlinear gravity theory with Lagrangian L = f(R), however, the energy current 
is not always divergence-free. From the field equation 

fiiRpv - -fg^u - V^V^/ij + g pu Of R = (8) 

where f R = df(R) / dR, one can find that the energy current is divergence- free only for the 
case with condition V p ^ v f R V p K u = 0. Moreover, one can not arrive at a true quasi-local 
energy from the integration (J4]) in general. Interestingly, for an FRW universe and static 
spherically symmetric solutions with constant scalar curvature, the quasi-local energy have 
been found [36(. For an FRW universe with the line element 



ds 2 



-dt 2 + a 2 (t)dp 2 + p 2 a 2 (t)d 2 n 2 



the generalized Misner-Sharp energy is 



E 



{pay 
12 



/ + 6Hf R - 6f R (H 2 + H) 



(9) 



Note that it is obviously not the special case of the mass-like function given in 21J , contrary 
to the case of GB gravity, see 
element 

ds 2 = -h(r)dt 2 + -^-dr 2 + r 2 d 2 n n _ 2 , (10) 
9{r) 

the energy for n = 4 is found as 



19]. For static spherically symmetric solutions with line 

-h(r)dt 2 + ~^dr 2 + r 2 d 2 tt n „ 2 , 
9{r) 



E= r - 
2 



f R -g} R + \r 2 U-f R R) 

o 



(11) 



where R, f R , and / are all needed to be constants. 
For the scalar-tensor gravity with Lagrangian 

L = F{$)R - \ (V$) 2 - V($) 



where F(<&) is an arbitrary positive continuous function of the scalar field $ and V($) is its 
potential, the equations of motion are 



FG^v - V M V,F + g^UF - 



1 



1. 



V At $V,$ - ^(-V A $ V A $ + V) 



8ttT, 



(12) 

□$ - V"($) + F'($)i? = 0. (13) 

In an FRW spacetime, the energy current is divergence-free since (V^V^F + | V^V^) V^K 1 
0. The generalized Misner-Sharp energy has been obtained as 



E 



(pa-Y 
12 



6 (fH 2 + HF} - - V 



(14) 



It has been proved that the unified first law, which was previously proposed in Einstein 
gravity, holds also for Gauss-Bonnet gravity, f(R) gravity of these two cases, and scalar- 
tensor gravity in an FRW spacetime, with the uniform 



V a E = A^ a + wv a v, 



(15) 



where A is the area of the sphere with radius r and V is its volume. W is called work density 
defined as W = -h ab T ab /2 and 



# a = T b a d b r + Wd a r 



(16) 



is the energy supply vector, with T ab being the projection of the n-dimensional energy- 
momentum tenor of matter in the normal direction of the (n — 2)-dimensional sphere. 



III. DYNAMICAL HORIZON ENTROPY 



In this section, we will follow Hayward's proposal 33j, |34j to define the dynamical horizon 
entropy directly by replacing the Killing vector with the Kodama vector in Noether charge 
method 



25 



261 ]. However, one should be careful that the usually used Wald entropy ex- 



pression (such as the one used in [28]) has invoked the property of Killing vector. In fact, 
just after Wald and Iyer presented their entropy expressions, Jacobson has pointed out 31] 
that for nonstationary case, there are three obvious candidate forms of the entropy, which 
are the full potential produced by the Noether charge construction associated with certain 
vector field, the potential eliminating the higher derivatives of the vector field via identities 
that would hold for the Killing vector, and the standard Wald entropy that is the potential 
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dropping the terms proportional to the vector field and the binormal of horizon is invoked 
to replace the first order derivative of vector. All of the three expressions yield identical 
results when pulled back to a bifurcate Killing horizon. But for a dynamical horizon, the 
three forms are different in general. Because the later two forms are derived involving some 
properties of Killing vector, we will consider the full Noether potential. In fact, if there is a 
fundamental criterion which can be imposed to define the entropy of a nonstationary black 
hole, after a deeper understanding of black hole entropy, the full Noether potential should 
be considered as the first candidate. 

For a generally covariant Lagrangian L, the variation of Lagrangian density n-form is 
described as 

where ipn) denote the field variables, E^ = is the equation of motion for ip^, and is 
a functional of the field variables and their variations. When we identify the variation with 
a general coordinate transformation induced by an arbitrary vector field <r, one can 

obtain 

Then the vector field = G M — <^L is divergence- free V^J^ = on shell, for which one 
can find an antisymmetric Noether potential satisfied with J M = V ' V Q^ V '. The full Noether 
potential Q^ v can be calculated in a straightforward manner for a given action, as shown 



m 



37. 



381 ] . For our aim, we will consider the Lagrangian involves no more than quadratic 



derivatives of metric g^ and the first order derivative of a scalar field $, so the action is 
given by 

1 = lib / dnx ^9 L ^u, R,ux p , V„$). (17) 
Variation of this action gives the equation of motion E^' and Q 13 as 

E { $ = M, u - Ig^L - X a ^R u)pPa - 2V P V A X (/ ^, 

E m = _ 



where 
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Regarding the variation as a coordinate transformation induced by an arbitrary vector q, we 
can obtain the Noether current 

J p = 2V a [X^Vfr - 2&V M X a ^] - 2^ - 2R avm X w ^ - ~u/V„$) . 

by manipulating the covariant derivatives and using the equation of motion. Since the 
general covariance of the Lagrangian implies 

(m? - 2R avpii X™^ - Vfrf = 

for an arbitrary vector 77^, the Noether potential of the current J p can be obtained as 

QT = -2X puXp V x <, p + ^ p V x X puXp . (18) 

Integrating the Noether potential over any closed spacelike surface B of codimension n — 2, 
the Noether charge is proportional to 

S= ~L! Q " udB ^ (i9) 



B 



where dB^ = ^e^-^d 71 2 y. When q p is a timelike Killing vector, the term proportional 
to q p of Q pv is absent in the integral, because q p vanishes on the Killing horizon. Then 



the entropy is reduced to the one used in 28]. For dynamical spacetimes, however, we 



will take q p as the Kodama vector, so this term must be preserved. Moreover, it should 
be noticed that for Einstein and Gauss-Bonnet gravity, while not for the f(R) gravity and 
scalar-tensor gravity, this term also disappears because \7\X puXp = in those cases. Hence 
one can expect more difference between the stationary and dynamical entropy for f(R) 
gravity and scalar-tensor gravity than for Einstein gravity and Gauss-Bonnet gravity. 

Actually, for Einstein gravity, Hayward has found that the entropy has same form as the 
stationary case. For Gauss-Bonnet gravity with Lagrangian (J5j), one has 

x ^Xp = g^ g \"\p] + 2a (g p[x g Mp] R + 2g" [x R Mp] - 2g p[x R Mp] + R puXp ) . (20) 



For simplicity, we consider only the case with n — 5. Using Eqs. (Ej) , (J3J) , (118]) and ( 1201) . we 

can evaluate the entropy (fl9|) . which is 

1 , 3o , 7r 2 r 3 q . , 

S = —A H —A = — h 6n ar. 21 

4 r 2 2 

One can find that it has the same form as the stationary case, and when a = it reduces 
to A/ A consistent with the entropy of Einstein gravity. 
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For f(R) gravity, we have 

Considering a four-dimensional spacetime, we obtain the dynamical horizon entropy 

S = nr f R - — . (22) 



For the scalar-tensor gravity with 



X^ xp = Fg 



lA. x n \ v \p\ 

y ) 



the entropy has the similar form as Eq. (|22 



S = 7l r 2 F-- r -^^. (23) 
4 r 



1UV 



It should be noticed that for the entropy expression given by the boost-invariant fields, the 



entropy of scalar-tensor gravity is generally different with the entropy of f(R) gravity 26 ]. 
But our result is reasonable since the f(R) gravity can be treated as a special scalar-tensor 
theory by introducing the scalar field <fi = R and potential V_j= <fif — f in the Brans-Dick 
theory, and choosing the Brans-Dick parameter u = (see 



for a review). 



The entropies ( 122]) and ( 123]) reduce to A/4 for Einstein gravity where fn = F = 1. It 
should be stressed that the stationary case has only the first term and the extra term may 
change our understanding about the nonequilibrium thermodynamics of f{R) gravity and 
scalar-tensor gravity in dynamical spacetimes, as we will show. 

Before doing that, we want to discuss whether or not the dynamical entropy satisfies the 
second law of thermodynamics. Hayward has shown that for Einstein gravity there is the 
second law 33] , which states that if the null energy condition holds on a future (respectively 
past) outer trapping horizon, or on a past (respectively future) inner trapping horizon, then 
the horizon entropy is non-decreasing (respectively non-increasing) along the horizon. We 
review the proof briefly. Denote the tangent vector to the horizon by £ = d/dX = j3d u — ad v 
and fix the orientations by 9 + = and (3 > on the horizon. Then = d9 + /d\ = 
(3d u 8 + — ad v 9 + yields dr/d\ = —otd v r = —f3r9^d u 9 + /2d v 9 + . Considering the null energy 
condition 

T uu > and T vv > 0, (24) 

and the Einstein field equation, one can know d u 9 + < 0. The signs of 6L and d v 9 + are given 
by the definition of future or past, outer or inner trapping horizons. Defining the directional 
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derivative along the tangent vector to the horizon 

5 X = D/dX = CV a , (25) 

one can obtain S X S = 2irrdr/dX > for future outer or past inner trapping horizons, and 
vice verse. 

Now we will give an alternative proof of the second law, which is important because it is 
applicable to modified gravity theories. Consider the unified first law ( TT51) projecting along 

cv a E = A* a e + wev a v. 

For Einstein gravity, it has been shown that this equation can be written as a Gibbs equation 

6 X E = ^6 X S + W5 X V. (26) 

Z7T 

Hence 5\S can be written as 

5 X S = 2irAV a £ a /K = -27vAe~ 4 '/3T uu d v r / ' k. 

Considering the null energy condition (12^|) and the surface gravity ([3]), one can immediately 
find 5 X S > for future outer or past inner trapping horizons, and vice verse. In the 
following sections, we will show that the Gibbs equation (126]) holds for Gauss-Bonnet gravity 
in any dynamical spacetime, and also holds for f(R) and scalar-tensor gravity in the FRW 
spacetime with slowly varying horizon. Thus our dynamical entropy expressions satisfy the 
second law in these cases. 



IV. EQUILIBRIUM THERMODYNAMICS OF GAUSS-BONNET GRAVITY IN 
DYNAMICAL SPACETIMES 

On any static horizon, the equilibrium thermodynamics of Gauss-Bonnet gravity has 
been constructed recently using Padmanabhan's identity [7|. For the dynamical case, one 
also has the equilibrium thermodynamics after projecting the unified first law along the 
horizon, where it was found that the entropy needs the same form as f T2"TT) . However, the 
thermodynamics is only restricted on FRW and Vaidya spacetimes js), In this section, 
we will study the generalization to any dynamical spacetime with spherical symmetry. 
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A. Hayward's identity 



We will first consider the method based on Hayward's unified first law. The components 
of tangent vector £ a can be determined by C, a d a d u r = up to a normalization of £ a , which 
is irrelevant for our aim. In double-null coordinates, we write it clearly 

C=(-d u d v r,d u d u r). (27) 

Using the field equation of Gauss-Bonnet gravity ([6]), the energy supply along the horizon 
is obtained as 

^^a^ [ T VV T )U (Tj U 4>m ~\~^mu ) ~\~ T-uvXiv ^"iuv ] 

3 71 " <f> 2 

T 1 iv 1 'iuv 1 'mu 37TQ!?")X) T^ uv Tf UU -\-0(t, u ), (^8) 

where 0(r, u ) denotes some terms proportional to r, u (or its higher orders). One can find 
that the variation of entropy (I2T!) along the horizon £ a V a S just equals to Eq. ( f28|) up to 
some different terms proportional to r, u and the temperature factor T = k/(2tt). So we can 
recast the unified first law along the horizon as a Gibbs equation 

5 X E = T5 X S + W5 X V. (29) 



B. Padmanabhan's identity 

Then we will extend Padmanabhan's identity to the dynamical spacetime. The identity 
has been constructed for Gauss-Bonnet gravity on static horizons, which reads 6] 

dE = TdS - PdV. (30) 



Although it is also a Gibbs equation similar to Eq. f[29|) . there are two key difference which 
should be clarified. First, the differentials d, which are different with 5\, are interpreted as 
dE = (dE/dr + ) dr + etc., where r + refers to the horizon radius. This means that we are 
considering two solutions to the gravitational field equations differing infinitesimally in the 
parameters such that horizons occur at two different radii r + and r + + dr + , instead of the 
case in Eq. ( 1291) where the variation is taken as the directional deriviative tangent to the 
horizon. So one must be careful that here all quantities should be evaluated on the horizon 
before manipulating the differentials. Second, in the static spacetime with line element f|T0|) . 

13 



the pressure in Eq. (130!) is interpreted as the radial pressure of matter on the horizon, 
given by P = T r r , which is seemly different with the work density W = — | (T/ + T r r ). In 
dynamical spacetimes with line element (JTJ), it is not known what is the radial pressure. We 
will use the work density as the desired pressure, because we will show its consistency with 
the (negative) radial pressure P on the horizon of static spacetime (fTUi) with n — 5. From 
the field equations of Gauss-Bonnet gravity, we have 

»< P? " ^ - \ Kr) *%$ mr) F + 4« - 4a 9(r) ] . 

Observing the regularity of the curvature scalars on the horizon, one must impose several 
restrictions on the metric functions 

g\r + ) = h'{r+), g{r + ) = h{r + ) = 0, (31) 

which leads to T r r = T\ on the horizon, i.e. P = —W on the horizon as expected. 
Thus our aim is to check the identity 

dE = TdS + WdV (32) 

on the dynamical horizon. The work density on the horizon can be derived as 

W = e*T uv = JL [ r + eVr, TO +4ae\ uv } r _ r . (33) 

Using the expression of entropy (1211) . one can easily obtain 

/ 3?r 2 r 2 \ 
dS = ( + 6n 2 a J dr + . (34) 



Using Eqs. (|3|h (|33|) , and (134|), we have 

TdS + WdV = -^dr + . 

AG 

The generalized Misner-Sharp energy ([7]) on the horizon is 

3vrri r „-, 
E = [1 + 2«r; 2 ] 



and its variation is 



dE = 3 ~^rdr + . (35) 



One can find that Eq. ( 132]) holds on the horizon actually. 

Obviously, the Gibbs equations (1291) and ( )32l) tell us that the equilibrium thermodynamics 
of Gauss-Bonnet gravity can be constructed from gravity field equations in any dynamical 
spacetime with spherical symmetry. 
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V. EQUILIBRIUM THERMODYNAMICS OF f(R) GRAVITY 



A. Hayward's identity 

Now we will study whether or not there is the Gibbs equation (1291) for f(R) gravity in 
an FRW spacetime. Using the coordinate transformation 

u = £* - p, v = t* + p (36) 

where is determined by dt*/dt = 1/a, the FRW spacetime can be described as the double- 
null form ([1]), with 

r — pa. (37) 



2 

In an FRW spacetime, the trapping horizon with 6 + = is located at 



P = ~, (38) 
a 



which coincides with the apparent horizon 



40J], that has the radius r + = 1/H. Using the 



transformation laws ( 1361) and ( 1371) . we can read the Kodama vector as 

K» = (I, -rH, 0,0), 

and the surface gravity on the horizon as 

k = -e+d tt d v r = -\{2H 2 + H) = -H(l + ^). (39) 

Consider the model with a = t q , which leads to k — —(q— that shows that the surface 

gravity is negative in general (except the radiation dominated case). So we are treating an 
inner trapping horizon, rather than outer trapping horizon. Using Hamilton- Jacobi method, 



it has been pointed out in the recent work 



41] that the temperature of the inner horizon 



is T = —k/ (2tt) preserving the positive temperature. So we need to add a minus before 
T in the Gibbs equations (129]) and (!32|) . which can be interpreted as that the energy may 
decrease when the entropy of inner trapping horizon increases. However, to be concise, we 
will still use T — k/ (2tt) for the inner horizon. 

The tangent vector along the horizon £ 6 (here index b = t,r) ( 1271) can be read as 

e = (l,~2 P H-^), (40) 
15 



up to a proportional factor. Using the field equation of f(R) gravity © and the tangent 
vector (j40p . the energy supply along the horizon can be obtained as 

b f R H(2H* + H) f R {2H 2 + H) } R (2H 2 + H) 
A ^ ~ ^4 + > ( 41 ) 



where we have used Eq. (1381) to restrict the result on the horizon. The entropy (1221) in the 
FRW spacetime can be expressed as 



4 W 2H 2 + H 

Multiplying the factor T = k/ (2tt) to the entropy variation along the horizon, we obtain 

_ f R H(2H 2 + H) _ f R (4H* + lQH 2 H + 3H 2 + 2HH) Jn_ 
X 2m AH%2H 2 + H) + 2fP' { ^ ] 

where we have used Eq. (138|) . It is interesting to find that Eq. (T4ll is same as Eq. (|42j) . 

provided that the horizon is varied so slowly that 

H < # 2 , F < if 3 . (43) 

So we have established the Gibbs equation (129]) for f(R) gravity. In the general case, 
however, Eq. (129]) does not hold. The difference between Eq. (f4T|) and Eq. ( l4"2j) can be 
given as 

Td H S = 5 A £ - W5 A V - ^5 A S 
= A* 6 £ 6 - ^hS 

271 

_ Hf R f R (6H 2 H + H 2 + HH) 
4# 4 4H 3 (2H 2 + H) 

B. Padmanabhan's identity 

Next we will check Padmanabhan's identity (|32|) . Before doing that, we will show that 
W is still consistent with — P for f(R) gravity in the static spacetime. This can be carried 
out by evaluating the difference between T\ and T r r 

(ti - r ;) = mimmmm^m^ + g{r)m = . (44) 



The last equality holds on the horizon since we have used Eq. ( |31j) . In an FRW spacetime, 
the work density on the apparent horizon can be get as 

W = ^ (/ - 6f R H 2 + 5Hf R - 4f R H + /r) . (45) 
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The key step is to consider the variations. For the apparent horizon of FRW spacetime, we 
notice that the differential d can be expressed as 

, d d dt d 

d = dr + — = dr + —^- = dr + — 1 - — . 

ar + a im a H(t) at 

Thus, we can replace the differential d with dt for the aim of checking Padmanabhan's 
identity ([32}. Using Eqs. ©, Q22J), (EH]), and (JSJ, one can obtain the right hand of Eq. 
(J32D as 

TdS + WdV ~ Td t S + Wd t V 

f H ( 5R2 + 3 #) /b ( 2jff4 + ^H 2 H + 4ij 2 + 



+ 



AH A 2H 2 

f R (2H 2 -H^ 
4# 4 



2# 3 (# 2 + 2ij 



(46) 



The left hand is 



iB~*B=-^ + ^ ^2^ + ^ (47) 

Comparing Eqs. ( I4"6"j) and (j4Tj) under the approximation (H3"l) . we have justified the Gibbs 
equation ( 1321) . In the case without the approximation, we have 

TdpS = dE- —dS - WdV 
2vr 

~ d t E - ^-d t S - Wd t V 

_ Hf R f R (6H 2 H + H 2 + HH) 
4# 4 AH 3 (2H 2 + H) 

Interestingly, one can find dpS ~ dnS, which suggests both of them have the same origin 
and Padmanabhan's approach is consistent with Hayward's one. 

We will further check Padmanabhan's identity ( 1321) in the static spacetime where the 
generalized Misner-Sharp energy is also found. We will not consider another identity based 
on Hayward's unified first law because the energy supply (I16p is vanishing on the horizon 
of the static spacetime (fTOl) . Note that its time component is vanishing obviously, and its 
radial component is \l/ r = T T r — |(T/ + T r r ) = for T\ = T£ on the horizon (pH]) . Now we 
will evaluate the right hand in Eq. (|30|) . The surface gravity is k = g'/2 and the entropy is 
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S = f R A/4. Using the field equation ([8]), we obtain 

t^mWW 1 2 f f , , f , r 2 gf' R h' , r 2 f R g'h> r 2 f R gh' 2 r 2 f R gh" 
TdS + WdV = -r f-rgf R + rf R g -^- + — _ + 

l if , f , , r 2 f R g'h' r 2 f R gh' 2 r 2 f R gh" 

\rlf + r + f R9 ' + r ^f^dr + , (48) 

where we have used f' R = in the second line with the mind that R, f, and f R are all 
constant, which is the requirement of quasi-local Misner-Sharp energy. Moreover, we have 
used Eq. (13 ip in third line. Reading the generalized Misner-Sharp energy ( TTTj) on the horizon 
and respecting that R, f, and f R are all constant, we can get the energy variation as 

dE = (±r 2 + f + l -f R - \rlf R R^j dr + . 

Substituting Ricci scalar on the horizon 

R _ 2(1 - g) 2g' 2gh' g'h' | gh' 2 gh" 

r 2 r rh 2h 2h 2 h 

2 V g'h" 

r 2 r h' ' 

where we have used Eq. (J3"T]) in the second equality, the variation dE can be recast as 

, P A,,, f , , r 2 + f R g'h" \ 
dE = ^-r + f + r + f R g + ^ j dr + , 

which is same as Eq. (j4"8"|) . Thus, we have shown that Padmanabhan's identity ( 1321) holds. 

VI. EQUILIBRIUM THERMODYNAMICS OF SCALAR-TENSOR GRAVITY 
A. Hayward's identity 

We will study the thermodynamics of scalar-tensor gravity in an FRW spacetime. Using 
the field equation (fl~2"|) and the tangent vector (l4Tfl) . the energy supply along the horizon can 
be obtained as 

b FH(2H 2 + H) F(2H 2 + H) F(2H 2 + H) &{2H 2 + H) 
A ^ = 2~m AIP + AH^ + 8^ • (49) 

Consider the second Friedmann equation 



1 

H =F 



-A^p + P )+ l -HF~ l -F--^ 2 
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Obviously, when H is small, the term HF/2 = H§dF/d§/2 ~ $ ~ F, is small in gen- 
eral. Then the last term in ( H9|) . which is a higher order small quantity, can be neglected. 
Moreover, we notice that 1/F takes roles as the effective Newton gravitational constant in 
the scalar-tensor theory. So $ 2 is very small indeed, since it is known that the experi- 
mental bounds acquire the effective Newton constant as an approximate constant |42|, i.e. 
F = <&dF/ d§ ~ $ is small. Read the entropy (123|) in the FRW spacetime as 

ff= 4( F __^), (50) 
4 V 2H 2 + H h V ; 

and its variation as 

rr, c _ FH{2H 2 + H) F(AH 4 + 16H 2 H + 3H 2 + 2HH) F 
6xb ~ 2~m AH\2H 2 + H) + ^ {b) 

One can find that the Gibbs equation (1291) holds under the approximation (H3|) . The differ- 
ence between Eq. fH9|) and Eq. (l5Tj) can be given as 

Td H S = 5 X E - W5 X V - ^5 X S 

= A* b e - ±S X S 

HF F(QH 2 H + H 2 + HH) ¥(2H 2 + H) 
~ ^H 1 + 4H 3 {2H 2 + H) + 8H 1 ' 

B. Padmanabhan's identity 

Next consider another Gibbs equation. Using Eqs. ( 1T2"|) and (|50p . one can obtain the 
right hand of Eq. fl32|) as 



TdS + WdV ~ T^S + Wd t V 

FH VH F (2H 4 + 13H 2 H + AH 2 + HH^j F (2H 2 - 
= ~2H 2 + Am 2R , ( 2H2 ~ ^ + 

The left hand is 

,„ orn FH f(h 2 -2H) p QVH + 3H& -2H$$ V 



(52) 



2H 2 2H 3 2H 2 2AH A 12H 3 

FH VH F[H 2 + 3H) J_ $ 2 (2H 2 + H) 



2H 2 AH 4 2H 3 2H 2 8H 4 

where we have invoked the equation of motion of scalar field (Tl3|) in the last equality with 
the mind that V = V'$. Comparing Eqs. fl52|) and fl53l) using the approximation f|43l) and 
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omitting the term with $ 2 , we have justified the Gibbs equation ( 132]) . In the case without 
the approximation, we have 

TdpS = dE- —dS - WdV 
2n 

~ d t E - ^-d t S - Wd t V 

271 

HF F(6H 2 H + H 2 + HH) &(2H 2 + H) 
' 4^4 + 4HH2IPTH) + 8H 1 ' 
which shows dpS ~ dnS. 



VII. CONCLUSION AND DISCUSSION 



In this paper, we have investigated the relationship between the gravitational field equa- 
tion and the thermodynamics on the dynamical trapping horizon with spherically symmetry. 
Following Hayward's proposal, we have obtained a general expression of Wald-Kodama en- 
tropy, by replacing the Killing vector as Kodama vector in the Noether charge method 
proposed by Wald. We have evaluated Wald-Kodama entropy in Gauss-Bonnet gravity, 
f(R) gravity and scalar-tensor gravity. It is shown that the dynamical horizon entropy 



of f{R) gravity and scalar-tensor gravity have different forms with t 



reir stationary cases, 



21 



24 



36 



43]. 



contrary to the assumption given in many references 0, 0, Q, Q, 

We have shown that the second law of thermodynamics is satisfied under the null en- 
ergy condition for the dynamical horizon entropy of arbitrary gravity theories, if Hayward's 
thermodynamic identity can be constructed using the entropy expression. We have built 
up this identity for Gauss-Bonnet gravity in any trapping horizon and for f(R) gravity and 
scalar-tensor gravity in the FRW spacetime with slowly varied horizon. Thus, the second 
law holds and we can argue that the dynamical entropy is preferred at least in these cases. 

Besides the Hayward's thermodynamic identity, we have found that Padmanabhan's ther- 
modynamic identity can be generalized to any trapping horizon for Gauss-Bonnet gravity. 
For f(R) gravity and scalar-tensor gravity, we have disclosed that the equilibrium identity 
can still be constructed in the FRW spacetime with slowly varied horizon. Hence our work 
provides an alterative method to absorb the nonequilibrium entropy production without 



introducing the mass-like function 19 



21 



43]. For a general FRW spacetime, it is still not 



known whether the nonequilibrium entropy production is necessary. We also study Pad- 
manabhan's thermodynamic identity for f(R) gravity in the static spherically symmetric 



20 



spacetime with constant scalar curvature. We are restricted in these cases since the gen- 
eralized Misner-Sharp energy was only found there. We have shown that Padmanabhan's 
equilibrium identity holds in the static spacetime. This supports the result given in 



28 



301 ]. which suggests that the nonequilibrium entropy production is dispensable at least for 



stationary cases. 
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